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Abstract. We report on the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state in two-dimensional rf-wave superconductors 
with magnetic field parallel to the superconducting planes. This state occurs at high magnetic field near the Pauli-Clogston 
limit and is a consequence of the competition between the pair condensation and Zeeman energy. We use the quasiclassical 
theory to self-consistently compute the spatially nonuniform order parameter. Our self-consistent calculations show that the 
FFLO state of a d-wave order parameter breaks translational symmetry along preferred directions. The orientation of the 
nodes in real space is pinned by the nodes of the basis function in momentum space. Here, we present results for the Knight 
shift and discuss the implications for recent nuclear magnetic resonance measurements on CeCoIns. 
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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
of spin-singlet superconductors is the compromise be- 
tween the pairing condensate, favoring anti-parallel spin 
alignment, and the Zeeman effect, favoring parallel spin 
alignment along the field (HBl- This compromise leads 
to a spatially inhomogeneous state of "normal" and "su- 
perconducting" regions, where the "normal" regions are 
defined by a spectrum of spin-polarized quasiparticles. 

The FFLO phase of d-wave superconductors is modi- 
fied by the anisotropy of the order parameter in momen- 
tum space compared to s-wave superconductors. The up- 
per critical transition line, B^2{T), has a kink at low tem- 
peratures, T* ^ 0.06 Tc, corresponding to the discontin- 
uous change in the modulation of the order parameter 
||2lSBl01- Recent calculations of the spatial modulation 
of the order parameter in 2D near Bc2 predicted that the 
energetically favored state at low T and high B forms a 
"square lattice" instead of the ID stripe order 

Here, we restrict our study to temperatures above 
this structural phase transition, T > T*, and address the 
quasiparticle response in the FFLO phase between the 
lower critical field Bd and the upper critical field Bc2- 
In addition, we assume that B is parallel to the super- 
conducting planes. In this geometry the magnetic field 
affects the superconducting condensate only through the 
Zeeman coupling of the quasiparticle spin to the field. 
Furthermore, we assume a cylindrical Fermi surface. 

Within the quasiclassical theory of superconductiv- 
ity we calculate self-consistently the order parameter 
A(R,p) and the quasiclassical Green's functions by solv- 
ing Eilenberger's equation in a constant magnetic field B. 
The Zeeman coupling of the quasiparticle spin with mag- 
netic field enters through /iB • C7, where (7, are Pauli spin 



matrices and ~ {g/2)pLB is the absolute value of the 
magnetic moment of a quasiparticle with negative charge 
e; pLB = \e\/2mc is the Bohr magneton. Note that the g- 
factor is a free material parameter in this calculation. 

From the solutions of Eilenberger's equation we can 
calculate measurable quantities like the free energy, 
quasiparticle local density of states 0] and local magne- 
tization M(R). Here, we consider spin-singlet order pa- 
rameters that factorize into A(R,p) = A(R) cos20, with 
a spatially dependent amplitude, A(R), and an angular 
dependent d^i_yi-wa\t basis function. 

The local magnetization is given by the paramagnetic 
response of the medium and the spin-vector component 
of the quasiclassical Matsubara Green's function | 



M(R) = IjxNf 
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with the normal-state density of states A^^ at the Fermi 
level. The normal-state susceptibility, Xn = 2/i^A^^, is 
defined by Mn = XnB. 

For comparison, we show in Fig. 1 the calculated tem- 
perature dependence of the magnetization in the uniform 
superconducting (USC) phase for three different values 
of B. Increasing B changes the T-dependence of the mag- 
netization from linear to quadratic with a residual zero- 
temperature value due to the field induced shift in the 
spin-split density of states of the gapless d-wavs super- 
conductor. This result is in agreement with scaling argu- 
ments by Yang and Sondhi lil(Xl . 

In Fig. 2, we show temperature scans of the minimum, 
average and maximum local magnetization for the stable 
FFLO phase, with spatial order-parameter modulation 
along the (110) direction, i.e., along the nodal direction 




FIGURE 1. Magnetization M of the spatially uniform super- 
conducting (USC) phase normalized by the normal-state mag- 
netization Mfti for different magnetic fields b = flB/lTzTe. 

of the gap function. The Knight shift K is proportional 
to the change in the local magnetic field at the nucleus, 
thus it is directly proportional to the local magnetization. 
Since K is weighted by the field distribution, the largest 
contribution comes from areas where the derivative of M 
vanishes, which are at the minimum and maximum loca- 
tions of M. The calculated bifurcation between minimum 
and maximum local magnetization seen in Fig. 2 is in 
qualitative agreement with measurements of the Knight 
shift on CeCoIns reported by Kakuyanagi et al. fllll . 

In Fig. 3, we show field scans of the local magnetiza- 
tion at r/r^ = 0.1 starting in the USC phase and into the 
FFLO (110) phase. It illustrates the nonlinear magnetic 
response of the quasiparticles due to an external field and 
the continuous second order transition at the lower crit- 
ical field Bel, which also is signaled by the appearance 
of a single domain wall. This fin ding clearly contradicts 
the claim by Yang and Sondhi llOll about a first order 
transition at B^i between the USC and FFLO phase. 

In addition, we calculated the spin-resolved local den- 
sity of states in the FFLO state [6] (not shown). We found 
that the characteristic Andreev bound states, due to the 
periodic sign change of the order parameter, are respon- 
sible for the excess spin polarization of quasiparticles at 
the domain walls seen in Figs. 2 and 3. Therefore, the 
Andreev bound states should be clearly visible features 
in scanning tunneling spectroscopy measurements. 
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FIGURE 2. Temperature scan of the normalized local mag- 
netization at fixed magnetic field b = 0.2 in FFLO phase for 
(110) orientation of spatial order-parameter modulation. Inset: 
Phase diagram with temperature and field scans across FFLO 
phase. 
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FIGURE 3. Field scan of the normalized local magnetization 
at fixed temperature T/r^ = 0.1 across USC and FFLO phases 
for (110) orientation of order-parameter modulation. 
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